Statics and dynamics of an incommensurate spin order in a geometrically frustrated 

antiferromagnet CdCr 2 04 
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Using elastic and inelastic neutron scattering we show that a cubic spinel, CdCr2 04, undergoes an 
elongation along the c-axis (c > a — b) at its spin-Peierls-like phase transition at Tjv = 7.8 K. The 
Neel phase (T < Tjv) has an incommesurate spin structure with a characteristic wave vector Qm 
= (0,5,1) with S ~ 0.09 and with spins lying on the ac-plane. This is in stark contrast to another 
well-known Cr-based spinel, ZnCr2 04, that undergoes a c-axis contraction and a commensurate spin 
order. The magnetic excitations of the incommensurate Neel state has a weak anisotropy gap of 0.6 
meV and it consists of at least three bands extending up to 5 meV. 

PACS numbers: 75.25.+Z, 75.30.Ds, 75.50.Ee 



When a transition metal oxide has degeneracy, orbital 
or magnetic, a novel phase transition can occur through 
the coupling of the relevant degrees of freedom to the lat- 
tice to lift the degeneracy. A well-known example is the 
Jahn- Teller lattice distortion that involves either doubly 
degenerate e g orbitals or triply degenerate ti g orbitals 
P, Q| ■ Geometrically frustrated magnets provide a fer- 
tile ground for similar novel phase transitions because 
the degeneracy can be macroscopic both for quantum 
and for classical spins Until now the most frustrat- 
ing system known is the one that consists of a network 
of corner-sharing tetrahedra, the pyrochlore lattice, with 
the simplest spin hamiltonian, Hi — J ^2 NN Si ■ Sj with 
isotropic antiferromagnetic nearest neighbor interactions 
only. Theoretically, it has been shown that in the ideal 
case, the spins alone cannot order even at zero temper- 
ature [2, Q . Experimentally the magnetic lattice can be 
realized in several materials such as in the pyrochlores 
AoB 2 7 @, spinels AB 2 4 and C15 Lave phases AB 2 
p. Among them, Cr-based spinels ACr 2 04 (A=Zn, Cd) 
realize the most frustrating lattice with the dominant 
antiferromagnetic nearest neighbor interactions due to 
the direct overlap of the ti g orbitals of the neighboring 
Cr 3+ (3d 3 ) ions HE] . Consequently, A&2O4 remains 
paramagnetic to temperatures far below the character- 
istic strength of the interactions between the spins, the 
Curie- Weiss temperature [6rrw| = 390 K and 88 K for 
A = Zn and Cd [II El Eg, respectively. Upon fur- 
ther cooling, however, the system undergoes a first or- 
der spin-Peierls-like phase transition Ih3 ] from a cubic 
paramagnet to a tetragonal Neel state at Tn = 12.5 K 
and 7.8 K for A = Zn 0] and Cd El El, respectively. 

Recently, ZnCr 2 04 has been studied extensively using 
neutron scattering techniques 0] . Its tetragonal dis- 
tortion involves a contraction along the c axis (c < a). 



Its Neel state has a rather complex commensurate spin 
structure. The spin structure has four different charac- 
teristic wave vectors, Qm, (5,5,0), (1,0,5), (5,5,5) and 
(0,0,1) 16]. Furthermore, the relative ratios of the neu- 
tron scattering intensities of these wave vectors vary de- 
pending on the subtle chemical conditions during sample 
preparation EH ■ This suggests that even in the tetrago- 
nal phase, ZnCr 2 04 is critically located close to several 
spin structures, that makes it difficult to understand the 
true nature of its ground state. The CdCr 2 04 compound 
in this class has not received much attention partly be- 
cause it was commonly believed that the same physics 
hold true as in ZnCr 2 04 E3, which we found is not the 
case. As the pyrochlore lattice possesses many ground 
states and multiple ways for lifting the ground state de- 
generacy, CdCr 2 04 provides a venue for understanding 
phase transitions in this seemingly complex class of ma- 
terials. 

In this paper, we report the results from elastic and in- 
elastic neutron scattering measurements on 11 CdCr 2 04 
(space group Fd3m, a = 8.58882 A for T — 10 K). 
Surprisingly, we find that it undergoes a phase transition 
that is qualitatively different in nature from the one ob- 
served in ZnCr 2 04. CdCr 2 04 elongates along the c-axis 
(c > a) and undergoes an incommensurate (IC) Neel or- 
der. The high Q-resolution data indicate that the incom- 
mensurate magnetic structure has a single characteristic 
wave vector of Qm = (0, 6, 1) with 5 ~ 0.09 perpendic- 
ular to the unique c-axis. We present two possible high 
symmetry spin structures that are consistent with the 
particular Qm- The interplay between the lattice dis- 
tortion and the IC spin structure is discussed along with 
dispersion of the spin wave excitations. 

Preliminary measurements were initially performed 
on the TAS2 thermal triple-axis spectrometer of Japan 
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FIG. 1: Color image of inelastic neutron scattering intensities 
from single crystals of CdCr2 04 at T = 15 K > Tjv with Tiuj 
— 0.6 meV. The data were taken at SPINS using eleven 2.1 
cm x 15 cm PG(002) analyzer blades with E; — 5meV. 



Atomic Energy Research Institute (JAERI), while more 
detailed studies ensued on the SPINS cold neutron triple- 
axis spectrometer of NIST Center for Neutron Research. 
A single crystal weighing ~ 100 mg was used for the 
elastic measurements, while three of these crystals were 
co-mounted within 1° mosaic for the inelastic measure- 
ments. The crystals were mounted in the (hkO) scattering 
plane, that allowed the investigations of three equivalent 
planes, (hkO), (Okl), and (hOl) planes due to the crystal- 
lographic domains. 

Fig. n shows the Q-dependence of the spin fluctua- 
tions in the spin liquid phase of CdCr204 measured above 
Tjv- The ring-shaped intensity around (2,2,0) is essen- 
tially identical to the one observed in ZnCr204, and it 
is due to the collective low energy excitations of antifcr- 
romagnetic hexagonal spin clusters in the pyrochlore lat- 
tice [l5|. This suggests that the cubic phase of CdCr204 
can be well represented by the Hamiltonian Tii as in 
ZnCr2C>4. J is estimated from the Curie- Weiss tempera- 
ture \@cw\ = 88.97 K (see Fig. 2(a)), to be J = -1.02 
meV. Although CdCr204 appears to have the same fun- 
damental spin degrees of freedom in the cubic phase as 
in ZnCr204, it exhibits strikingly different behaviors in 
the tetragonal phase below Tv = 7.8 K. It undergoes 
an elongation along the c axis (Fig. |2|b) and (d)) and 
the magnetic long range order has an incommensurate 
characteristic wave vector (see Fig. Etc)). 

The positions of the magnetic Bragg reflections found 
in the scattering plane are shown in Fig. Et a )- The possi- 
ble characteristic incommensurate wave vectors are Qm 
= (0,(5,1) or (0,1,5) or (1,(5,0) with 5 = 0.0894(3). (h and 
k are interchangeable.) However, since c > a, the three 



FIG. 2: (a) Bulk susceptibility, x, as a function of T. The 
transition is of first-order and has a small but finite hysteresis 
over 0.1 degrees. The inset shows a linear fit to 1/x- (b) T- 
dependence of the normalized integrated neutron scattering 
intensity at the magnetic peaks and that of the lattice strains, 
£ a, c = -£r£y~- The data were obtained from fitting the data in 
(c) and (d) to gaussians. (c), (d) Elastic neutron scattering 
data (c) through magnetic (1,5,0) IC points and (d) through 
a nuclear (400) Bragg reflection below and above Tjv. 



possible Qjvr's would produce the magnetic Bragg reflec- 
tions at slightly different positions in the scattering plane. 
We performed elastic scans with a high Q-resolution over 
a set of three IC positions to distinguish between the dif- 
ferent scenarios. Fig. 3(b) shows the results. The black 
solid, blue dotted and red dashed arrows correspond to 
the expected peak positions for Qm = (0, <5, 1), (1, (5, 0), 
and (0,1,(5), respectively. All three of the experimental 
peak positions are consistent with Qm = (0,(5,1), which 
indicates that the incommensurability occurs either along 
the a- or the 6-axis, perpendicular to the elongated c- 
axis. Furthermore, we have performed polarized neutron 
diffraction at TAS1 of JAERI (the details and results 
of the experiment will be reported elsewhere |l7j). and 
found that the spins are lying on the plane that is perpen- 
dicular to the incommensurability direction. For conve- 
nience, we chose Qm = (0,(5,1) with the spins lying in the 
ac-plane. The resulting breakouts of the magnetic Bragg 
reflections into three crystallographic and two magnetic 
domains are shown as different symbols in Fig. 3 (a). 

Following group theory arguments 18], the magnetic 
lattice with the characteristic wave vector of Qm consists 
of four independent sublattices as represented by spheres 
with different colors in Fig. Et c ) and (d). Each sublat- 
tice connects every third nearest neighboring Cr 3+ ions 
that are separated by the symmetrically equivalent dis- 
tance of (5,5,0). These are the second nearest neighbors 
along the chains represented by lines in Fig. Etc) an d 
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FIG. 3: (a) The observed pattern of the magnetic reflec- 
tions on the scattering plane. The quadruplets consist of two 
doublets at Qjv±(0,5,1) and Qjv±(5,0,1), respectively. The 
central peaks of the triplets are tails of the doublets that are 
present above and below the scattering plane. The colored 
symbols distinguish the crystallographic domains, (b) Elastic 
Q-scans centered at the three peaks under a grey shadow in 
(a). (c),(d) two possible spin structures deduced from Qm 
= (0,(5,1). Spins are rotating on the ac plane. The double 
and single lines represent the nearest neighbor bonds in the 
basal ab plane with J a b and in the out-of-plane with J c , re- 
spectively. In (c) the NN spin orientations along J c are close 
to perpendicular while in (d) they are close to collinear. 



(d). Within the sublattice, spins are aligned according 
to Sj = S e 27 "QA''( r J~ r °) and therefore the neighbour 
spins rotate by 2a, 7r, and 2a + 7r in the ab-, ac- and fee- 
plane, respectively, where 2a = 5 ■ it = 16.2°. In order to 
construct the relative orientation between the different 
sublattices, the magnetic interactions in the tetragonal 
phase of CdCr204 are considered. Analysis of a series of 
chromium oxides indicates that dJ/dr « 40 meV/A [l9j |. 
This implies that the tetragonal distortion (c > a = b) 
yields stronger AFM interactions in the basal plane with 
J a b = —1-19 meV (double lines in Fig. 3 (c) and (d)) 
and weaker AFM interactions among all other spin pairs 
with J c = —0.95 meV (single lines). For each chain in the 
basal plane, the nearest neighboring spins that belong to 
two different sublattices would favor a phase difference of 
7r ± a to minimize the exchange energy due to the strong 
AFM Jab- The out-of-plane J c is frustrating that allows 
stacking of the < 110 > chains along the c axis. Fig. Etc) 
and (d) show two highly symmetric spin structures that 
have the same mean-field exchange energy; in one, the 
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FIG. 4: Constant-?itj ((a)-(c)) and constant-Q ((d)-(e)) scans 
of the spin wave excitations in CdCr2 04. The solid lines are 
guide to eyes. 



chains are stacked almost orthogonally along the c axis 
(Fig. Pic)) and the other, the chains are stacked almost 
collinearly (Fig. E2d)). 

To understand the elementary excitations of the IC 
Neel state, constant-?kj and contant-Q scans were per- 
formed (Fig. The constant-Q scan at an IC zone 
center (Fig. Eld)) shows that there is a gap, A, of 0.65 
meV along with at least two additional excitation peaks 
around 2.3 and 4.7 meV. At Q values away from the zone 
center, the lowest energy peak shifts considerably in en- 
ergy while the two higher energy peaks shift only slightly 
(Fig. E](e)-(f)). Fig- El( a )-( c ) show the dispersiveness 
at low energies. Additional scans were performed to map 
out the dispersion relation of the magnetic fluctuations 
along the (l,k,0) direction (the arrow from Qy = 0.5 to 
2.5 in Fig. EI a ))? summarized in Fig. EI e )- 

What kind of spin Hamiltonian would select the ob- 
served IC Neel state as the ground state of CdCr 2 04 
in the tetragonal phase? The exchange anisotropy of 
the nearest neighboring (NN) couplings, J c ^ J a0l alone 
would favor a commensurate spin structure rather than 
the IC one, which means that additional perturbations 
to Tii are present. However, A ~ 0.6 meV is small 
compared to the entire energy band width of the dis- 
persion, ~ 5 meV, that suggests that the additional per- 
turbations must be small. Let us first discuss the ef- 
fect of the exchange anisotropy of the NN couplings, 
J c 7^ J a b, on the dispersion of the two model spin struc- 
tures shown in Fig. 3 (c) and (d). In the ideal isotropic 
case, J c = Jab, there is a four-fold (six- fold) degenerate 
flat zero energy mode and another four-fold (two-fold) 
degenerate dispersive mode that extends to higher en- 
ergies in the orthogonal (collinear) model (Fig. Et a ))- 
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FIG. 5: (a), (b), (c), (d) The spin wave calculations with 
the orthogonal (solid lines) and the collinear (dashed lines) 
commensurate models with different ratios of J c /J a b- In (a), 
(b), and (c) the energy is in unit of J a i,S while in (d) the 
energy is calculated in meV with J a t, = —1.19 meV and J c = 
—0.95 meV. Since the overlap between the (l,fe,0) and (0,fc,l) 
domains is considered, the k = 1 point corresponds to the 
AFM zone center in the former case while it corresponds to 
the zone boundary for the latter, (e) The dispersion relations 
determined from constant-Q and contant-7kj scans. The line 
is guide to eyes. The data were taken along the direction in 
the Q space shown as an arrow in Fig. |3f a). The open and 
closed symbols are the data obtained with focusing and flat 
analyzers, respectively. 



As the anisotropy is introduced, the zero-energy mode 
becomes dispersive. In the collinear structure, the six- 
fold zero-energy mode splits into two-fold and four-fold 
weakly dispersive modes, whereas it remains four-fold de- 
generate for the orthogonal structure j2jj. The disper- 
sions are not affected by the exchange anisotropy at high 
energies as much as they do at low energies. Indeed, other 
anisotropy terms such as the single ion anisotropy and the 
Dzyaloshinsky-Moriya (DM) interactions also modify the 
dispersion at low energies considerably but not at high 
energies. Fig. Efd) shows the dispersion relations ob- 
tained using the experimentally determined J a b = —1.19 
meV and J c = —0.95 meV. This shows that the orthogo- 
nal spin structure fits the dispersion at high energies bet- 
ter than the collinear one. This means that the Neel state 
of CdCr204 is likely to be the incommensurate spin struc- 
ture with the orthogonal stacking shown in Fig. |3fc). 



The actual incommensurability seems to be caused by 
other perturbative terms, such as further nearest inter- 
actions and/or DM interactions. We also considered the 
exchange interactions between third nearest neighbor in- 
teractions, and extensively examined the phase space as 
a function of J. We found a region where an incommen- 
surate spin structure can be selected as a ground state 
but the incommensurability was along the < 110 > direc- 
tion, as recently observed in LiCuVC>4 |2l|, rather than 
the observed < 010 > direction. This suggests that the 
tetragonal distortion involves distortions of the oxygen 
octahedra that lowers the crystal symmetry. 

In summary, we identified a spin-lattice coupling mech- 
anism that lifts the magnetic frustration in CdC^CU 
and that is distinctly different from the one observed in 
ZnC^O-j. The tetragonal distortion involves an elon- 
gation along the c-axis and the Neel state has a he- 
lical spin structure with the single characteristic wave 
vector of Qm = (0,(5,1). Our identification of the spin 
structure and dynamics of the low temperature phase of 
CdCr2C>4 should provide a unique test to theoretical at- 
tempts to explain the spin-Peierls-like phase transitions 
in the Heisenberg pyrochlore antifcrromagnets. 
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